NACA TN 3811 


\ 


03 

NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 3811 


CHARTS ADAPTED FROM VAN DRIEST’S TURBULENT FLAT -PLATE 

THEORY FOR DETERMINING VALUES OF TURBULENT AERODYNAMIC 

FRICTION AND HEAT -TRANSFER COEFFICIENTS 

By Dorothy B. Lee and Maxime A. Faget 

Langley Aeronautical Laboratory 
Langley Field, Va. 



Washington LIBRARY CM 

October 1956 

OCT 26 1956 


LANGLEY AERONAUTICAL LABGRATOf- 
, . LIBRARY NACA 

" r Y FIELD, VIRGINIA 







NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 38ll 


CHARTS ADAPTED FROM VAN DRIEST' S TURBULENT FLAT-PLATE 
THEORY FOR DETERMINING VALUES OF TURBULENT AERODYNAMIC 
FRICTION AND HEAT-TRANSFER COEFFICIENTS 
By Dorothy B. Lee and Maxime A. Faget 


SUMMARY 


A modified method of Van Driest' s flat-plate theory for turbulent 
boundary layer has been found to simplify the calculation of local skin- 
friction coefficients which, in turn, have made it possible to obtain 
through Reynolds analogy theoretical turbulent heat-transfer coeffi- 
cients in the form of Stanton number. A general formula is given and 
charts are presented from which the modified method can be solved for 
Mach numbers 1.0 to 12.0, temperature ratios 0.2 to 6.0, and Reynolds 

numbers 0.2 x 10^ to 200 x 10^. 

INTRODUCTION 


Van Driest' s theory has been found to be the most convenient of the 
turbulent theories as a means of estimating skin temperatures for the 
many flight trajectories at the Langley Pilotless Aircraft Research 
Division for various Mach numbers, temperature ratios, and Reynolds 
numbers; the theory has agreed well with experimental data. However, 

Van Driest 's equation cannot be solved directly and a rapid means of 
obtaining theoretical values was desired. Inasmuch as the presentation 
of Van Driest 's theory in the form of plots in sufficiently small incre- 
ments for accurate and easy interpolation would require too many plots 
to be practical, a modification was devised which, with one equation and 
three charts, gives values of skin-friction coefficients in the range of 
Mach numbers 1.0 to 12.0, temperature ratios 0.2 to 6.0, and Reynolds num- 
bers 0.2 x 10^ to 200 x 10^. This modified method along with the nec- 
essary charts is presented herein. 


SYMBOLS 


c f 


2t 

local skin-friction coefficient, — E 

' pV2 


specific heat of air at constant pressure, BTU/lb-°F 
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g 

gravitational force, 32.2 ft/sec 2 

h 

aerodynamic heat-transfer coefficient, BTU/sec-ft 2 -°F 

k 

correction factor 

M 

Mach number 

%t 

Stanton number , — ^ — 

gpVc p 

R 

Reynolds number , 

T 

temperature, °R 

V 

velocity, ft/sec 

X 

characteristic length, ft 

Z 

slope of straight-line approximation of c^ against R 
(fig. 1) 

7 

ratio of specific heats, 1.4 

P 

density of air, slugs/cu ft 

T 

shear stress, slugs/ft-sec 2 

P 

absolute viscosity of air, slugs/ft-sec 

(JO 

constant in power law for viscosity, namely, p w = pJ^i/ 


Subscripts : 


w 

wall 

00 

free- stream condition 

S 

straight line 
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DISCUSSION 


Method for Turhulent Flow 


Van Driest' s flat-plate theory for turbulent boundary layer (ref. l) 
in the form assuming the Von Karman similarity law for mixing length 
(ref. 2) is expressed in the equation for the local skin friction as 


0.242 

clri“l 

A - B 

2 A + gln -l 


B 

w 

A(= f a/ 2 fe:V /2 ' 



(~ 

I 2 + 1 

V f >°°) \T m ) 


\2A) J 

[_\ 2 A/ 

1 


- 11/2 


O.la + log 10 R„c fj00 - 0 ) log 10 


( 1 ) 


where 


7 - 1 M 2 
2 2 

hr = — £ 


1 + M„ 2 

B = d 1 

Tw/Tco 


0) = 0.76 


The value of Cp x obtained from equation (l) may be used to deter- 
mine the heat-transfer coefficient in the form of Stanton number according 
to reference 3‘* 


N St ~ °*^ c f 


( 2 ) 


Some values of cp were obtained from equation (l) for a range of 

m 5 

— from 0.2 to 6.0, Mach numbers from 1.0 to 12.0. and R from 0.2x10 

rp ' ' 

x co 

r 

to 200 x 10° with a punched-tape digital computer making as many as 
11 iterations . The curves of Cp against R on log plots for given 

Mach numbers and temperature ratios were found to closely approximate a 
straight line faired through the values for Reynolds numbers of 1 x 10 ^ 

and 4-0 X 106, as can be seen from the typical plots in figures l(a) 
and l(b) which are for Mach numbers of 1 and 5, respectively, and for 
T 

— of 0.2, 2.0, and 6.0. This straight line can be expressed 

Iqo 


k 
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' f > S ' R Z 


(5) 


and. since s is equal to at R = 1 x 10 °, 


A = ( C A 6 l 1 x 10 V 

V f /R=lxlO° 


Thus, 


{*) 

'r=ixio 6 

/ R .) 

r 

li x io6j 


00 


In order to evaluate Z it is necessary to substitute another set of 
values into equation (4), such as R = 40 x 10 ^ and Cf at R = 40 x 10 ^ 

/ c \ . = ( Cf )R=lxio 6 

\ f /R=A0xl0 6 fkO x 10 6 \ Z 
\ix io 6 j 

and 

log ( Cf )R=lxlQ 6 ., 

z = r f )R=40xl0 6 ( 5 ) 

log 4-0 

/* 

It can be seen that, with a chart of c^» at R = 1 x 10 and a plot of 
Z as a function of M and values of c f s can be obtained easily 

•^oo y 

T 

from equation (4) for any conditions of M, and R. 


NACA TN 3811 


5 


Values of Z for a range of Mach numbers from 1.0 to 12.0 at tem- 
perature ratios of 0.2 to 6.0 have been computed from equation ( 5 ) by- 

using values of c f at R = 1 x 10^ and R = 40 x 10^ obtained from 

the solution on the digital computer of Van Driest' s equation. (See 
eq. (l).) The values of Cf are listed in table I. The values of Z 

are shown in figure 2 (a) plotted against Mach number for several temper- 

T 

ature ratios. A cross plot of Z against — for several Mach numbers 

^•00 

is shown in figure 2 (b) which can conveniently be used to aid in inter- 
polation. Accurate values of c f calculated by the digital computer 

for a value of R of 1 x 10^ for use in equation (4) are shown in fig- 
ure 3(a) plotted against Mach number and are cross-plotted in 3(b) 

against . 


Because Van Driest' s curves of c f against R depart from the 

straight-line expression of equation (4), a correction factor k depend- 
ent on R is necessary for a higher degree of accuracy. The local skin- 
friction coefficient can then be written 

c f = kc f ^ s ( 6 ) 


Values of k have been determined by comparing c f s with solutions 
of Van Driest 's equation for Mach numbers 1, 3 , 6 , 9 , and 12, at temper- 
ature ratios of 0.2, 2 . 0 , and 6.0, and at Reynolds numbers of 0.2 x 10^, 

4 x 10 , 10 x 10 y and 200 x 10 . Figure 4 gives the correction factor 
k against R for the three temperature ratios and can be used for all 
Mach numbers from 1.0 to 12.0 since k proved to be almost invariant 
with Mach number. 

The present method consists in solving equation (4) for c^ by 
use of figures 2 and 3 and then determining Cf from equation ( 6 ) by 
use of figure 4. This procedure yields values of c^ well within 
±1.0 percent of Van Driest 's theory. If c f s is not corrected for 
k (fig. 4), the value obtained will be within ± 3.0 percent of the theory 
for R between 0.5 x 10^ and 100 x 10^. 
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Laminar Flow 

As a matter of convenience a chart is included for the determination 
of Ngt for laminar-flow conditions. Figure 5 is a plot of Ng^^R^, 

T 

against M for various values of obtained from reference 4. 

00 


ACCURACY 


£> 6 

The basic values of c^ at R = 1 x 10 and R = 40 x 10 
(figs. 3(a) and 3(b) and table I) , which were computed from Van Driest 's 
equation (eq. l) by the punched-tape digital computer with several iter- 
ations (approximately ll), are correct to four significant figures. 

The curves of Z (figs. 2(a) and 2(b)) are faired smooth from 
point to point and when used in equation (4) will give essentially exact 

agreement at R = 1 X 10^ and R = 40 x 10^ . 

The maximum error in k due to neglecting the influence of M is 
±0.75 percent. 


CONCLUDING REMARKS 


A modified method of Van Driest 's flat-plate theory for turbulent 
boundary layer has been found to simplify the calculation of local skin- 
friction coefficients which, in turn, have made it possible to obtain 
through Reynolds analogy theoretical turbulent heat-transfer coefficients 
in the form of Stanton number. A general formula is given and charts 
are presented from which the modified method can be solved for Mach num- 
bers 1.0 to 12.0, temperature ratios 0.2 to 6.0, and Reynolds numbers 

0.2 x 10^ to 200 x 10 6 . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June 20, 1956. 


NACA TN 3811 


7 


REFERENCES 


1. Van Driest, E. R.: 
Jour. Aero. Sci., 


Turbulent Boundary Layer in Compressible Fluids, 
vol. 18 , no. 3, Mar. 1951, pp. 14-5-160, 216. 


2. Van Driest, E. R.: The Turbulent Boundary Layer With Variable Prandtl 

Number. Rep. No. AL-1911, North American Aviation, Inc., Apr. 2, 1954 . 

3- Rubesin, Morris W.: A Modified Reynolds Analogy for the Compressible 

Turbulent Boundary Layer on a Flat Plate. NACA TN 2917, 1953. 

4. Van Driest, E. R.: Investigation of Laminar Boundary Layer in Com- 

pressible Fluids Using the Crocco Method. NACA TN 2597, 1952. 


8 


NACA TN 3811 


TABLE I 

VALUES OF SKIN-FRICTION COEFFICIENT 


M 

T 

ta = 0.2 

Too 

hi = 0.6 

Too 

T 

hi = 1.0 

Too 

^ = 2.0 
Too 

T ' 

=£ = 3.0 ; 

Too ' 


Ftf 3 

II 

ON 

b 

c f at R = 1 X 10 6 

0-5 

0.005100 

0.004170 

0.003675 

0.002990 

0.002610 

0.002350 

0.002005 

1.0 

.004904 

.004070 

.003603 

.002948 

.002578 

.002328 

.002000 

2.0 

.004290 

.003709 

.003345 

.002800 

.002475 

.002251 

.001950 

3-0 

.003596 

.003258 

.003005 

.002595 

.002328 

.002137 

.001875 

k.o 

.002978 

.002812 

.002652 

.002357 

.002153 

.002000 

.001779 

5.0 

.002475 

.002417 

.002325 

.002126 

.001974 

.001854 

.001675 

6.0 

.002078 

.002084 

.002036 

.001907 

.001798 

.001708 

.001566 

7-0 

.001764 

.001809 

.001788 

.001711 

.001635 

.001568 

.001460 

8.0 

.001516 

.001581 

.001582 

.001539 

.001489 

.001441 

.001357 

9.0 

.001317 

.001394 

.001407 

.001389 

.001357 

.001323 

.001262 

10.0 

.001156 

.001237 

.001258 

.001259 

.001241 

.001218 

.001173 

12.0 

.000913 

.000995 

.001025 

.001047 

.001046 

.001039 

.001017 


c f at R = 40 X 10 6 


.5 

.003033 

.002383 

.002048 

.001607 

.001370 

.001210 

.001000 

1.0 

.002915 

.002321 

.002007 

.001583 

.001350 

.001195 

.000997 

2.0 

.002537 

.002107 

.001857 

.001499 

.001293 

.001154 

.000970 

3.0 

.002113 

.001838 

.001660 

.001382 

.001211 

.001091 

.000930 

4.0 

.001736 

.001577 

.001456 

.001250 

.001115 

.001017 

.000879 

5.0 

.001433 

.001346 

.001267 

.001120 

.001016 

.000938 

.000824 

6.0 

.001193 

.001152 

.001102 

.000999 

.000921 

.000860 

.000767 

7.0 

.001006 

.000993 

.000962 

.OOO89O 

.000833 

.000785 

.000711 

8.0 

.000859 

.000862 

.000845 

.000796 

.000754 

.000717 

.000658 

9.0 

.000742 

.000754 

.000746 

.000714 

.000683 

.000655 

.000608 

10.0 

.000646 

.000666 

.000663 

.000643 

.000620 

.000599 

.000562 

12.0 

.000505 

.000530 

.000533 

.000528 

.000517 

.000505 

.000482 



(a) Mach number of 1. 

Figure 1.- Variation of local skin-friction coefficient with Reynolds 

number . 
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(b) Mach number of 5* 
Figure 1.- Concluded. 
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(a) Variation with Mach number for several temperature ratios. 

Figure 2.- Values of Z as a function of Mach number and temperature 

ratio. 
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(t,) Variation with temperature ratio for several Mach numbers 

Figure 2.- Concluded. 





(a) Variation with Mach number for several temperature ratios. 

Figure 3 .- Local skin-friction coefficient for a Reynolds number of 
1 X 10^ as a function of Mach number and temperature ratio. 
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(b) Variation with temperature ratio for several Mach numbers. 

Figure 3- “ Concluded. 
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Figure 4.- Correction factor for all Mach numbers. 
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Figure 5.- Local heat-transfer coefficient for laminar boundary layer 
of a compressible fluid flowing along a flat plate. 
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